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Landau levels and Riemann zeros 
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The number N(E) of complex zeros of the Riemann zeta function with positive imaginary part less 
than E is the sum of a 'smooth' function N(E) and a 'fluctuation'. Berry and Keating have shown 
that the asymptotic expansion of N(E) counts states of positive energy less than E in a 'regularized' 
semi- classical model with classical Hamiltonian H — xp. For a different regularization, Connes has 
shown that it counts states 'missing' from a continuum. Here we show how the 'absorption spectrum' 
model of Connes emerges as the lowest Landau level limit of a specific quantum mechanical model 
for a charged particle on a planar surface in an electric potential and uniform magnetic field. We 
suggest a role for the higher Landau levels in the fluctuation part of N(E). 



At the begining of the 20th century Polya and Hilbert 
conjectured that the imaginary part of the complex zeros 
of the Riemann zeta function £(s) are the eigenvalues of 
a self-adjoint operator H. The existence of such an oper- 
ator implies the celebrated Riemann hypothesis that all 
complex zeros lie on the 'critical line' IZe(s) = 1/2. Sup- 
port for this 'spectral' approach did not emerge until the 
1950s, when Selberg found a 'trace' formula relating the 
eigenvalues of the Laplacian on a compact hyperbolic sur- 
face to its geodesies; there is a strong resemblance of this 
formula to the Riemann- Weil 'explicit' formula relating 
the Riemann zeros to the prime numbers. Further sup- 
port came in the 1970s and 80s from the works of Mont- 
gomery and Odlyzko: assuming the Riemann hypothesis, 
they showed that the local statistics of the Riemann zeros 
is described by the Gaussian Unitary Ensemble (GUE) 
of Random Matrix Theory (see e.g. (H). Inspired by 
this, and by analogies between trace formulas of number 
theory and the Gutzwiller formula for chaotic dynamical 
systems, Berry conjectured the existence of a classical 
chaotic Hamiltonian system with isolated periodic orbits, 
with periods related to the prime numbers, such that the 
spectrum of the quantum theory gives the complex Rie- 
mann zeros . The GUE statistics of the zeros requires 
the Hamiltonian to break time-reversal invariance. 



Berry's conjecture received support in 1999 from the 
work of Connes [§| , and Berry and Keating 0] , on a semi- 
classical model for a particle in one dimension with classi- 
cal Hamiltonian H — xp. This Hamiltonian breaks time- 
reversal invariance since (x,p) — > (x,—p) =>• H — > —H. 
Classical orbits in this model are unbounded hyperpolas 
in phase space but they may be converted into closed or- 
bits by an identification at the boundaries of a restricted 
region of the phase plane that defines some 'regulariza- 
tion' of the model. This regularization also makes finite 
the number of states with energy less than E, and this 
number is related to the number N(E) of complex zeros 
of the Riemann zeta function with positive imaginary 
part less than E. The Riemann-van Mangoldt 'counting' 



formula states that 

N(E) = N{E) + ilm log C Q + iE^j (1) 

where N(E) is a 'smoothed' counting function that gives 
the average number of zeros, which is corrected by the 
remaining 'fluctuation' term. The smooth term can be 
written as 



N(E) = 



0(E) 



1. 



where (see e.g. 



H ( E)=in, io,,r [ - -log*-, 
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which is the phase of the zeta function on the critical line. 
The smooth term has the following asymptotic expansion 
for large E: 
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while the fluctuation term is of order log E. In the Berry- 
Keating model this asymptotic expression for N(E) is 
recovered from a semi-classical count of states with pos- 
itive energy less than E; the correct 0(1) term is ob- 
tained by taking into account the Maslov phase implied 
by the identifications that close the hyperbolic orbits 0] • 
In the Connes model, which has a different 'regulariza- 
tion', the same asymptotic expression is found as a semi- 
classical count of states missing from a continuous spec- 
trum. Thus, the Berry-Keating interpretation of N(E) 
is as an 'emission' spectrum while the Connes interpre- 
tation is as an 'absorption' spectrum. 

The first aim of this letter is to show how these semi- 
classical results can be recovered by viewing the H = xp 
model as a lowest Landau level (LLL) limit of a quan- 
tum mechanical model for a charged particle on the xy- 
plane in a constant uniform perpendicular magnetic field 
of strength B, and an electric potential ip = — Xxy, where 
A is a constant. Let fi be the particle's mass and — e its 
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charge. In the Landau gauge A = Bxdy for the vector 
potential 1-form A, the particle's Lagrangian is 

L = t;{ ±2 + V 2 ) - — V x ~ eXxy. (5) 
2 c 

This Lagrangian is invariant under (x,y) — > (— x, — y), 
which implies a double degeneracy of states of energy E. 
It is also quadratic, with two normal modes: the standard 
'cyclotronic' mode of the A = model, and a 'hyperbolic' 
mode. The respective angular frequencies are 

eB eB 
u> c = — cosh??, Uh — i — sinh$, (6) 

flC [IC 

where sinh(2$) = (2A/ic 2 )/ (eB 2 ); for simplicity, we shall 
consider the limiting case uj c ^> \u>h\, for which 
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In the quantum theory, the energy in the cyclotronic 
mode is (n + l/2)huj c for integer n labelling the Landau 
level. At sufficiently low energy only the lowest (n = 0) 
level is relevant and the physics is effectively described 
by the Lagrangian (see e.g. 0) 
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where £ is the 'magnetic' length. This LLL model is the 
(unregularized) one-dimensional H = xp model. From 
the Landau model perspective, each quantum state is as- 
sociated to a quantum of magnetic flux $ = 27r he/ e that 
occupies an area 2ir£ 2 in the xy-plane. In the LLL limit, 
this implies the standard semi-classical quantization in 
which one quantum state is associated with an area 2irh 
in phase space. 

The total number of states in each Landau level, in 
particular the lowest one, is infinite. To make the number 
of states in each level finite, we shall restrict the particle 
to the square 



\x\<L, 



\y\<L. 



(9) 



This implies for the LLL model a restriction on the phase 
space equivalent to that proposed by Connes for the 
H = xp model. It also implies that there is a maximum 
magnitude for the classical energy: \E\ < L 2 /£ 2 in units 
for which hwh = 1 ■ which fixes the value of A in terms of 
e and B by ([7]). For large L/l, the semi-classical estimate 
for the total number of quantum states of any energy al- 
lowed by this bound is L 2 /2tt£ 2 in each of the four quad- 
rants in the xy-plane separated by the x and y axes. The 
classical trajectory for a particle of energy E < E max is 
the hyperbola xy = E£ 2 , so the available classical phase 
space is the region in the square that lies between the 
two branches of the hyperbola. The number N SC (E) of 
semi-classical quantum states with energy less than E is 



the area of this region divided by 2ir£ 2 , but because of 
the double-degeneracy due to the (x, y) — > (—x — y) sym- 
metry we count only those states with (x > 0,y > 0), 
which yields 



E 



The first term diverges as L — > oo and may be inter- 
preted as a regularization of the continuum of states in 
the lowest Landau level for a particle free to move on 
the entire xy-plane. The finite correction, which has a 
magnitude equal to the asymptotic approximation of ^ 
to the smooth part of the counting formula, is negative, 
which led Connes to interpret it as representing spectral 
lines missing from the continuum. 

Another aim of this letter is to use the Landau model 
perspective to arrive at a proper quantum-mechanical un- 
derstanding of the regularized H — xp model of Connes; 
we focus on that model (details of which may be found 
in a recent book 0]) because we have not seen how to 
achieve the same objective for the Berry-Keating model. 
The Hamiltonian operator corresponding to the starting 
Lagrangian ([5]) is 



H 



2/2 



[fx + (Py + px) 2 } + eXxy , 



(11) 



where p x = —ihd x and p y = —ihd y . There is a unitary 
transformation in which this Hamiltonian becomes the 
sum of a Hamiltonian H c for the cyclotronic mode and 
a Hamiltonian for the hyperbolic mode. In units for 
which h = £ = 1 , one finds that 



H h = l -^(QP + PQ), (12) 



where P = —id/dQ and p = —id/dq. Note the Weyl 
ordering in Hh, which follows from the starting Hamil- 
tonian. The unitary transformation that achieves this 
decomposition simplifies considerably in the limit that 
u) c |wfc|; it corresponds to the classical canonical trans- 
formation 

q = X+p y , P = Px, Q = ~Py, P = V + Px- (13) 

In the limit that u c ^> \u>h\, any low energy eigenstate 
is the product of the ground state of H c , which is the 
gaussian ip — exp(—q 2 /2£ 2 ), with an eigenstate of H h , 
which we may choose to be either even or odd under 
Q — * — Q. These eigenstates of definite 'Q-parity' are 
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l<3l" 2 



The corresponding wavefunctions of the original Hamil- 
tonian are then given by 

^±(x,y) = C [ dQe-W 2 e -^ 2 / 2£2 $|(Q) . (15) 



FIG. 1: Absolute value of the wave function ipg(x,y) for 
E — 10 in the interval —10 < x,y < 10. The picture on 
the left gives a 3D representation while that on the right is a 
density plot. 



where C is a normalization constant. This yields 

'1 iE 1 (x-iy) 2 
J + T' 2' 2l 2 , 

lM*,fO - C* s (x- iy )e - + V- 



^(x,y) = C+e ^M(^ + 



(16) 
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where M(a,b,z) = iFi(a,b, z) is a confluent hypergeo- 
metric function. A representative plot of the even wave- 
function is shown in Fig [I] 

A classical trajectory that starts at (x, L) ends at 
(£,x), so that identification of these points creates pe- 
riodic orbits. This identification also means that the 
square in the xy-plane in which the particle is allowed 
to move is topologically a sphere, as shown in Fig. [2l al- 
though classical orbits lie entirely in one of four segments 
of this sphere. 

For these periodic orbits to emerge from the quan- 
tum theory in a semi-classical limit, we must identify 
the wavefunctions at these points, up to a phase. To 
determine the phase we use the asymptotic formulas 



M(a,b,z)~{ %o 



EM. c z v a-b Re ( z ) > o 
^ (-*)- Re(z) < 



(17) 



to deduce that 

iPe(L,x) ~ e" 
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L 2 \ -l/i-iE/2 



(18) 
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(19) 



FIG. 2: Left: icy-plane with the classical trajectories xy = 
El 2 . Right: Sphere that results from the identifications of the 
edges. Now the open trajectories become closed trajectories 
in one of four segments, enclosing one of the equatorial points 
A,B,C,D. 



for the odd functions. Observe that the 'even' wave- 
functions are symmetric, and the 'odd' wavefunctions 
antisymmetric, under the symmetry operation (x, y) — > 
(—x, — y). Counting both symmetric and anti-symmetric 
functions corresponds, semi-classically, to counting states 
for xy > 0, rather than just (x > 0, y > 0), so we could 
define the quantum model by considering only the sym- 
metric wavefunctions, which would effectively mean that 
we identify (x,y) with (— x, — y). The above formulas 
suggest that we should impose the boundary condition 

r E {x,L) = e *W^+-(,-i)/4 ^ e{L)X)) v = ±. (20) 

Applied to the 'even' (r; = +) energy cigcnfunctions, this 
leads to the asymptotic condition 



r(i + *f) (l 



-iE 



r(i-if) \2P 

which is equivalent to 

2i8(E) ( J±_\ 

\2ttPJ 



= 1 



(21) 
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where 0(E) is the function of fl3|), related to N(E) by j2]). 
This condition implies that 
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N(E) = N E 



(23) 



where N E is an integer that we may identify, asymptoti- 
cally, with the number of states with energy less than E. 
If the first term on the left hand side is interpreted as the 
rcgularization of the infinite number of states in a con- 
tinuum, then we see that N(E) has the Connes interpre- 
tation as a count of states missing from this continuum. 
This analysis can be repeated for the odd wave function 
tp E (x,y), in which case the function 9(E) is replaced by 
the phase of the odd Dirichlet L-functions. 
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Finally, we conclude with some speculations on how 
the fluctuation part of the Riemann counting formula 
might arise in the Landau model approach. In the con- 
text of the Berry-Keating model, the region of phase- 
space allowed by the regularization may be adjusted to 
'fluctuate' in such a way that the number of states of 
positive energy less than E is precisely N(E) [8j. How- 
ever, the Berry-Keating model does not seem to arise as 
the LLL limit of a model for a particle on a plane, and 
the 'fluctuating boundary' idea of |8[ does not work for 
the Connes model. A natural possibility for the Connes 
model, now viewed as a lowest Landau limit, is to sup- 
pose that the fluctuation term in the Riemann counting 
formula is related to the higher Landau levels. An imme- 
diate drawback of this idea is that the full Landau model 
is two-dimensional whereas we need a one-dimensional 
model. We thus seek some one-dimensional limit of the 
Landau model that generalizes the standard LLL limit. 
To this end, we return to the Hamiltonian as the sum of 
H c and Hf ll as given in (|12p and rescale, for convenience, 
to arrive at the Hamiltonian 



finds the semiclassical formula 



H = l(f + <,>) + 



PQ 



(24) 



where 7 = u> c /\u>h\- Introducing the standard annihi- 
lation operator a = (l/v / 2)(g + ip) and corresponding 
creation operator, such that [a, a 
the operators 



1, we observe that 



A = aQ l 



A ] = a'Q- 



(25) 



have a similar commutation relation but commute with 
the Hamiltonian: 



[A,rf]=l, [A,H]=0=[A\H] 



(26) 



Eigenstates of A are coherent states with complex eigen- 
values z. States with z = are those of the lowest Landau 
level, annihilated by a, so we may modify this limit by 
considering eigenstates of A with non-zero eigenvalue. At 
the classical level, the equation A = z implies that the 
cyclotronic motion is governed by the hyperbolic one. 
Taking z real for simplicity, one finds that 



q = V2zcos(7log<5), 



P 



-V2zsin( 7 logQ) . (27) 



Proceeding semi-classically, we consider the area of phase 
space enclosed by a closed classical trajectory, i.e. its 
action. The total action receives contributions from both 
modes: 



A 



PdQ+ pdq 



(28) 



In the limit, 7 >> 1 the first term on the RHS of this 
equation gives the estimate (|10j). Hence in the same limit 
and allowing the particle to gyrate according to (|27p , one 
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The first term, which diverges as L — > 00, may be in- 
terpreted as a count of states that become a continuum 
in this limit. Note the O(log-E) correction to the count- 
ing of 'missing' states, which is exactly the order of the 
fluctuation term in the Riemann counting formula. 

In summary, we have shown that the semi-classical 
H = xp model, conjectured to be the semi-classical limit 
of a quantum mechanical model for the complex zeros 
of the Riemann zeta function, can be understood as the 
lowest Landau level limit of a fully quantum mechanical 
model for a particle on a plane in the presence of electric 
and magnetic fields. As it stands, a counting of states in 
the model does not yield the full 'counting' formula for 
the Riemann zeros because the 'fluctuation term' is miss- 
ing, but we have provided semi-classical evidence that 
this will arise from a consideration of the higher Landau 
levels. Apart from providing a new tool in the spectral 
approach to the Riemann hypothesis, there is also the 
possibility that it will allow a laboratory construction of 
a system for which the physics is described by the 'Rie- 
mann Hamiltonian', the existence of which would prove 
the Riemann hypothesis. 
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